Probing the topological exciton condensate via Coulomb drag 
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The onset of exciton condensation in a topological insulator thin film was recently predicted. We 
calculate the critical temperature for this transition, taking into account screening effects. Further- 
more, we show that the proximity to this transition can be probed by measuring the Coulomb drag 
resistivity between the surfaces of the thin film as a function of temperature. This resistivity shows 
an upturn upon approaching the exciton-condensed state. 



Introduction. — Recently, there has been great experi- 
mental and theoretical interest in a new class of materials 
called topological insulators (TIs) pQ. These materials 
combine, for a particular doping level, an insulating be- 
havior in the bulk with topologically protected conduct- 
ing surface states that are described by a two-dimensional 
(2D) Dirac-Weyl Hamiltonian. When the surfaces of a TI 
film are independently doped or gated, electrons can be 
induced in one layer and holes in the other. In 2009, Ser- 
adjch et al. [2] predicted that such a system can support 
a so-called topological exciton condensate. Electrons in 
one layer combine with holes in the other layer to form 
excitons, which condense for low enough temperatures. 
One of the interesting features of this state is that it 
supports vortices in the order parameter: these vortices 
carry a fractional charge ±e/2 0. 

In a drag experiment with a two-layer system a current 
is applied through one of the layers (the drive or active 
layer) and an induced voltage is measured in the other 
layer (the passive layer). The drag resistivity prj is de- 
fined as the ratio between the electric field in the passive 
layer and the current density in the active layer. In the 
absence of an ordered state, the low-temperature behav- 
ior of the drag resistivity is /?d oc T 2 [or T 2 log(T) for a 
2D system with short-range interactions] . This quadratic 
temperature dependence was first observed in 1991 by 
Gramila et al. ]3j and is a hall mark of Fermi-liquid be- 
havior. A review of drag effects in two-layer systems is 
given in Ref. [4] . A departure from the standard pu oc T 2 
temperature dependence indicates the occurrence of non- 
Fermi-liquid behavior. The low-temperature dependence 
of pd for semiconductor electron-hole bilayers was stud- 
ied experimentally in numerous recent works [SJ |B] . On 
the theoretical side, we note the calculation by Hu [7], 
which predicted an enhancement of Coulomb drag for or- 
dinary electron-hole bilayers at temperatures just above 
the excitonic instability, and the analysis of Ref. [5], 




FIG. 1: The drag resistivity pu [in units of pr> = (uoU) 2 h/e 2 , 
with vo the density of states at the Fermi level and U the 
interlayer interaction strength] calculated from Eq. Q as a 
function of the reduced temperature T/T c . With this choice 
of units pd/pd depends on two parameters - the chemical 
potential p and the critical temperature T c - which implicitly 
determine the interaction strength U in Eq. |l]). The curves 
correspond to the case of equal chemical potentials in the 
two layers (p — 0.04 eV). From top to bottom the critical 
temperatures are Tc/Tp = (1/150,1/200,1/250). The vertical 
dashed lines are located at T = T c (left) and at the "upturn" 
temperature scale (see text) T = T u (right). 

which predicted a sharp increase of the drag resistivity 
as the temperature is lowered below the critical temper- 
ature of condensation. 

In this Letter we show that the drag resistivity shows 
a precursor of the topological exciton condensed phase, 
via an upturn upon approaching the critical tempera- 
ture, as shown in Fig. [T] Close to the critical tempera- 
ture we find that on the basis of a Boltzmann analysis 
pu cx log(T — T c ). This upturn could therefore be used 
to determine the proximity to the phase transition. We 
note that the system considered here resembles double- 
layer graphene (DLG) - a system of two graphene layers 
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FIG. 2: The Bethe-Salpeter equation for the effective inter- 
action V c ff - Eq. |T]). The thin wavy line is Vb and the ar- 
rows are the electron propagators in the top layer (top ar- 
rows) and the bottom layer (bottom arrows). Furthermore, 
K = fei — k$ = — is the center-of-mass momentum of the 
electron-hole pair and q = fci — k^ = ks — h.2 is the interlayer 
momentum transfer, which is responsible for Coulomb drag. 
Finally, p is an integration variable. To make contact with 
the drag resistivity in Eq. (|3| replace k = ki and k' — fc2. 



separated by a dielectric barrier. The main differences 
are the number of degenerate electron species and the 
dielectric constant [HI EH- The results presented here 
are therefore also qualitatively applicable to DLG, a sys- 
tem for which exciton condensation has been predicted 
too [21 [HI Q2] . Coulomb drag in DLG has been inten- 
sively studied theoretically |13j . This literature, however, 
refers to drag in which both layers are either electron- or 
hole-doped so that exciton condensation does not occur. 

Coulomb drag and exciton condensation. — We con- 
sider a TI thin film whose top layer is hole-doped and 
whose bottom layer is electron-doped. We apply the 
"closed-band approximation" [SJ [TTJ [T2], i.e., we con- 
sider only the upper Dirac cone for the electron layer 
and the lower Dirac cone for the hole layer. This is jus- 
tified since the temperatures we consider are sufficiently 
small and the screening lengths sufficiently large, so that 
the far-lying bands can in first instance be considered 
inert. Then, the dispersions are well approximated by 
eb(fe) = v\k\ — /ib and et(fe) = — v|fc| + m for the bot- 
tom and top layer, respectively (h — 1 throughout this 
Letter). Here v ss 5x 10 5 m/s is the Fermi velocity appro- 
priate for the surface states of a typical TI thin film [T3] . 
These formulas are valid up to a cutoff £ = 0.2 eV, which 
is the typical distance between the Dirac point and the 
bulk bands in a TI. We also introduce the mean Fermi 
energy, /i = (/ib + Mt)/2, and half the chemical potential 
imbalance, h = (/ib — Mt)/2- 

Our theoretical treatment is based on the Boltzmann 
equation, whose main merit is to give a physically trans- 
parent picture of the scattering processes that control 
the momentum transfer rate between the layers. The cru- 
cial quantity in the Boltzmann approach is the scattering 
amplitude V e s (ki, k%, ki), which is diagrammatically 
presented on the left-hand side in Fig. [2j where fci and 
&2 are incoming momenta, one from each layer, and k% 
and ki outgoing momenta. 



Under ordinary circumstances this interaction is well 
approximated by the single screened interaction line 
Vb(fci — ki), which is shown as a wavy line in Fig. [2] 
But this simple picture breaks down in the vicinity of 
the exciton pairing transition, where it becomes neces- 
sary to include the effect of pairing fluctuations, dia- 
grammatically represented by the infinite series of lad- 
der diagrams shown in the right-hand side Fig. [2j This 
series of diagrams diverges at the critical temperature 
T c when the center-of-mass momentum of the electron- 
hole pair K = k\ — k$ and the corresponding energy 
fl = eb(fei) — ^t(ks) tend to zero. The series can be 
calculated analytically if the momentum dependence of 
the screened interaction Vo is neglected, which is a rea- 
sonable approximation in most cases and sufficient for 
a qualitative understanding of the effect of pairing on 
Coulomb drag. One then obtains an effective interaction 
that depends only on the center-of-mass momentum K 
and energy as 

V(k!, k 2 , k 3 , fc 4 ) ~ V e s(K, Q) = — ^ , (1) 

where U is the momentum-independent contact interac- 
tion strength representing the screened interaction and 



1 V n(e t (k + K)) - n(e h (k)) 
A ^ e b (fe) - e t (k + K) - Q. - i0+ 



(2) 



is a pairing susceptibility, A being the 2D electron sys- 
tem area. The temperature T enters the above ex- 
pression through the Fermi-Dirac distribution, n(E) = 
1/[1 + exp(/3.E)], where /3 = {k^T^ 1 . It is easy to check 
that, in the symmetric case (h — 0) S(0, 0) diverges 
logarithmically at zero temperature, and therefore the 
denominator of Eq. ([!]) must vanish at some non-zero 
temperature T c no matter how small the electron-hole 
attraction. This is the mean-field pairing transition tem- 
perature. We see that, as the temperature is decreased 
towards T c , pairing fluctuations first lead to an enhance- 
ment of V e ff for small K and 17 and ultimately to a pole 
at fi = and K = as T = T c . The occurrence of 
this pole is the cause for the upturn of the Coulomb drag 
resistivity close to T c . 

Inserting Eq. ([T]) in the Boltzmann collision integral 
and performing standard manipulations, we finally arrive 
at the following expression for the drag resistivity py>: 



Pd = 



2(2tt) 



dKdfl 



\V cS (K,fl)\ 2 
sinh 2 (/?fi/2) 



dkdk'^m [b(k; K, fi)] 3m [b[k'; K, Q)} 



x [v t (k' +K)- v t (k + K)] ■ K(fc') - v h {k)] , (3) 

where v t ^(k) = Ve t (b)(fe) are the group velocities, n is 
the single layer carrier density, and 

n(e t (k + K))-n(e h (k)) 
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(4) 
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is the summand of Eq. Since v t (k) and Vb(k) are 
oppositely directed po is positive, as expected for carriers 
of opposite polarities. 

Drag Resistivity. — To determine the drag resistivity 
we need to evaluate Eq. ^ numerically. However, the 
qualitative behavior near T c can be obtained analytically 
as follows (in what follows we first consider the balanced 
case h = 0). First, we notice that the denominator of 
Eq. fl) can be expanded, for small K and f2 as follows: 



1 - UE(K,0) 



a(T) + a(T)Uv (f3vK) 2 
iUv Q pn/A , 



(5) 



where a(T) cx T — T c measures the distance from the 
critical point, vq — n/(2irv 2 ) is the density-of-states at 
the Fermi energy, and a(T) is a (dimensionless) posi- 
tive ultraviolet-convergent quantity. Further, for K > 
we have [v t (k' + K) - v t (k + K)] ■ [v h {k') - v h (k)} = 
2v 2 [cos((j>) — 1], where <fi is the angle between k and k' . 
Since this quantity vanishes only in a set of zero measure 
with respect to the two 2D integrals over k and k' in 
Eq. ([3]), we can approximate the scalar product with a 
momentum-independent quantity of the order of v 2 . The 
integrals over k and k' can now be carried out analyt- 
ically: using that J dk Qm[b(k; 0, fl)] — Sm[S(0,f2)] = 
— fo/30/4 for pVt < 1, we obtain 



dKdfl- 
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{v U) 2 K 



(6) 



which is immediately seen to diverge logarithmically 
when T — > T c (i.e., for a — > 0). 

In Fig. [l] we show the numerically evaluated drag re- 
sistivity pd as a function of temperature. We observe 
that for temperatures much larger than T c , the drag re- 
sistivity increases quadratically as po oc T 2 15j . The 
logarithmic divergence of pd oc log(T — T c ) is clearly 
visible. From top to bottom the curves correspond to 
T C /T F = (1/150,1/200,1/250). Following Ref. [6], we 
define T u as the temperature at which the upturn of 
prj starts, as shown in Fig. [T] As a rule of thumb, 
T u « (3/2)T e . 

Critical temperature. In order to assess whether 
the upturn in the drag resistivity can be observed exper- 
imentally, we need to calculate the (mean-field [16]) crit- 
ical temperature as a function of interlayer distance and 
carrier density, taking into account screening effects |17j . 
We employ a separable approximation to the momentum- 
dependent interlayer interaction Vb(fei — k^) [TH]. The 
advantage of this approach over using a contact inter- 
action is that it captures the decrease of the interaction 
strength Vq with increasing transferred momenta. Fur- 
thermore, no ultraviolet cutoff on which the result for T c 
would depend is needed. 



The screened interlayer Coulomb interaction is given 
by Vo(q,w) = V t h{q)/e{q,Lo), where V t h(q) is the bare in- 
teraction between one electron on the top surface and one 
on the bottom surface, q — \k — k'\ and u) = eb(k) — eb(k') 
are momentum and energy transfers, respectively, and 
e(q,uj) is the dielectric screening function of the carriers, 
which we approximate in the random phase approxima- 
tion. Both Vtb(<z) and e(q,u>) depend on the dielectric 
constants of the TI and of the environment surrounding 
the thin film. The details of the construction are given 
in Ref. [T§]. 

We consider a TI film with vacuum above the top sur- 
face and a typical substrate material (e.g. SiG^) below 
the bottom surface: the appropriate background dielec- 
tric constants are etop = 1, £ti = 100, ebottom = 4. 
For the balanced case with h = we obtain the phase 
diagram in Fig. [3] From right to left the lines corre- 
spond to carrier densities n — (0.25, 0.5, 1) x 10 10 cm~ 2 . 
Our results show that the critical temperatures are within 
reach using existing cryogenic techniques for solid-state 
systems. Furthermore, also the required surface carrier 
densities have been reached previously pQ. In the inset 
of Fig. [3] we show the effect of a non-zero density imbal- 
ance. The solid curves represent the behavior of T c /p 
for constant kpd. From top to bottom the curves corre- 
spond to k F d = (0.05, 0.075, 0.1). We note that T c /n de- 
creases strongly with increasing density imbalance. The 
dashed line is the line along which the quartic term of the 
Landau-Ginzburg expansion of the free energy in powers 
of the order parameter, evaluated at T c , vanishes. Thus, 
for a fixed h-pd, the transition becomes first order as h 
grows beyond the intersection of the dashed line and the 
curve corresponding to that kpd value. The first-order 
region below the dashed line cannot be accurately de- 
scribed by our present normal-state formalism. The de- 
creasing behavior of T c with h indicates that the excitonic 
transition is strongly suppressed by the presence of any 
density imbalance. Thus, it is important in experimen- 
tal realizations to have accurate control over the electron 
and hole densities in the two layers. 



Discussion and conclusions. — Our theory predicts a di- 
vergence of pd as T approaches T c , but does not include 
the possible suppression of this divergence by critical fluc- 
tuations in the immediate vicinity of T c . We expect that 
inclusion of these fluctuations leads to a down turn of 
Pd for temperatures very close to the critical tempera- 
ture [2U], before increasing again upon further decrease 
of the temperature, as predicted in Ref. [8]. 

We note that if the layers are brought too close to- 
gether, the wave functions in the top layer and the bot- 
tom layer begin to overlap and interlayer tunneling be- 
comes important. This limits our theory to a minimal 
interlayer distance d > 10 nm. Actually, also conduc- 
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FIG. 3: Critical temperature T c (in K) versus interlayer dis- 
tance d (in nm) in the absence of density imbalance (h = 0). 
From right to left the lines correspond to carrier densities 
n = (0.25,0.5, 1) x 10 10 cm" 2 . Inset: The critical lines T c /p 
versus h/p. From top to bottom the lines correspond to 
fcpd = (0.05,0.075,0.1). The critical lines terminate on the 
dashed line, i.e., the locus of the points where the fourth order 
term of the Ginzburg-Landau expansion of the free energy in 
powers of the order parameter vanishes. 

tion over the side layers of the TI film could be impor- 
tant. One expects that for samples of large enough area, 
these can be neglected. If not, a transport-prohibiting 
gap could be induced in the side layers by interfacing 
those with a ferromagnet. 

In closing, we note that the upturn of p^> is analo- 
gous to the upturn of the spin-drag resistivity recently 
predicted in a cold gas of fermionic atoms near a ferro- 
magnetic transition [21 . In particular, spin-fluctuation 
contributions to the effective interaction in the ferromag- 
netic case play the same role as pairing fluctuations in 
the present case. The reason why in Ref. [ST] the spin 
drag resistivity was predicted to remain finite at T c is 
that the transverse spin-fluctuation channel was treated 
as part of the familiar (Hartree) charge and longitudi- 
nal spin-density fluctuation channels. When transverse 
spin fluctuations in the Fock channel are also included, a 
logarithmic divergence is obtained, just as in the present 
case. We also note that our results, while qualitatively 
similar to those obtained in Ref. [7] for massive electrons 
and holes in a conventional electron-hole semiconductor 
bilayer, differ, at the theoretical level, by the inclusion 
of an additional term in the collision integral. This term 
is easily missed by using the Kubo formula. The nature 
of the additional term is similar to the vertex corrections 
that are responsible for replacing the momentum lifetime 
by the transport lifetime in the classical Drude formula 
for the resistivity. 
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